#we start with the Y matrix
wine.sensory <- matrix(0,8,5)

#row.1 <- cbind(6,5,3,4,2)
wine.sensory[1,] <- cbind(6,4,3,4,2)
wine.sensory[2,] <- cbind(5,3,2,3,2)
wine.sensory[3,] <- cbind(7,6,3,2,2)
wine.sensory[4,] <- cbind(5,2,1,4,4)
wine.sensory[5,] <- cbind(5,3,2,4,6)
wine.sensory[6,] <- cbind(7,1,4,5,8)
wine.sensory[7,] <- cbind(6,2,3,6,9)
wine.sensory[8,] <- cbind(5,1,4,4,9)
rownames(wine.sensory) <- c('Merlot','Cabernet','Shiraz','Pinot','Chardonnay','Sauvignon','Riesling','Gewurtz')
colnames(wine.sensory) <- c('Fruity','Spicy','Sweet','Acidic','Hedonic')

#Let's now make the X matrix.
wine.chemical <- matrix(0,8,3)
wine.chemical[,1] <- rbind(533,514,516,437,434,600,770,670)
wine.chemical[,2] <- rbind(138,139,143,135,133,133,123,124)
wine.chemical[,3] <- rbind(275,241,220,300,261,317,215,251)
rownames(wine.chemical) <- rownames(wine.sensory)
colnames(wine.chemical) <- c('Acidity','Alcohol','Sugar')





